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The obhque parameters S, T, and U and their higher-order extensions (V, W, and X) are observ- 
ables that combine electroweak precision data to quantify deviation from the Standard Model. These 
parameters were calculated at one loop in the basis-independent CP-violating Two-Higgs Doublet 
Model (2HDM). The scalar parameter space of the 2HDM was randomly sampled within limits 
imposed by unitarity and found to produce values of the oblique parameters within experimental 
bounds, with the exception of T. The experimental limits on T were used to predict information 
about the mass of the charged Higgs boson and the diflerence in mass between the charged Higgs 
boson and the heaviest neutral Higgs boson. In particular, it was found that the 2HDM predicts 
-600 GeV < ra^± — mz <100 GeV, with values of mu± > 250 GeV being preferred. The mass 
scale of the new physics (Mnp) produced by random sampling was consistently fairly high, with 
the average of the scalar masses falling between 400 and 800 GeV for Y2 = nr^, although the model 
can be tuned to produce a light neutral Higgs mass (~ 120 GeV). Hence, the values produced for V, 
W, and X fell well within .01 of zero, confirming the robustness of the linear expansion approxima- 
tion. Taking the CP-conserving limit of the model was found to not significantly affect the values 
generated for the oblique parameters. 



1. INTRODUCTION 

In supersymmetric models, the scalar sector of the theory features multiple Higgs bosons. In the minimal 
SUSY scenario, which requires two Higgs doublets, the phenomenology features three neutral Higgs bosons 
and one charged Higgs boson pair. However, if multiple Higgs particles are discovered in the absence of 
evidence of squarks or other SUSY partners, it will not be immediately apparent what symmetries gov- 
ern these phenomena. Because a generic Two-Higgs-Doublet-Model (2HDM) produces unrealistically large 
flavor-changing neutral currents, naturalness arguments would lead one to expect that some of symmetry is 
present in the Higgs sector, but the symmetry need not be associated with SUSY. In this paper we examine 
a generic 2HDM in a formalism that does not impose symmetries on the scalar sector, as in refs. pTl and 
Absent such symmetries, all physical observables must be basis- independent (invariant with respect to 
transformation in Higgs- flavor space). 

In the case that the scale of the new physics is not too much larger than the Z boson mass, the 2HDM 
produces non-zero shifts of the oblique parameters S, T, and U as well as the higher-order parameters V, 
W, and X. In this paper we calculate the six parameters in the 2HDM and randomly sample the 2HDM 
parameter space to compare the predictions of the model to the experimental results. The goal of this work 
is to present basis-independent formulae for the extended parameters (V, W, and X) and to answer the 
following questions: How do the values generated in the 2HDM for the oblique parameters compare with 
the experimental values? Does the 2HDM predict significant shifts in the higher-order parameters? (ie, is 
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any additional information provided by going beyond S, T, and U?) What do these resuhs for the obhque 
parameters reveal about the phenomenology (neutral and charged scalar masses) of the 2HDM? 

We also examine factors that may affect the validity of our results, such as the differences between the 
various definitions of the oblique parameters appearing in the literature, the effect of our choice of I2 (an 
unconstrained parameter in the scalar potential) on the size of the scalar masses, the errors inherent in 
using experimental limits from the PDG that use a fixed reference mass, and the extent to which taking the 
CP-conserving limit changes the numerical results. 



1.1. The Basis-Independent Two-Higgs Doublet Model 

The theory of a basis-independent CP-violating 2HDM has been developed in refs. [2] and [3]. What 
follows is a brief summary of the elements of the theory that will be relevant for the calculations in this 
paper. 

The most generic scalar potential for a model of two Higgs doublets is 

V = i;s$i<f b + ^z,-,,M<^,){<i>l<^>d) , (1.1) 

where the Higgs doublet fields are written as ^a{x) = ($^(x) , <i>°(a;)), with a = 1, 2; complex conjugation 
converts a barred index to unbarred and vice versa. The fields $a can be redefined by an arbitary U(2) 
transformation which rotates them to a different basis in Higgs flavor space. The coefficients Y^i and Z^j^^j^ 
in eq. (1.1) are likewise basis-dependent. Any measurable parameters of the 2HDM must be invariant with 
respect to such a transformation. Invariants can be constructed by summing over all Higgs flavor indices, 
pairing unbarred and barred indices, as is done in the covariant equation |1.1[ 
The vacuum expectation values of the Higgs doublets are 

i^^) = ^A?A^ (1-2) 



\/2 V ^« 

where v = 246 GeV and Va is a vector of unit norm. One can deflne new doublet fields as follows: 

= {H+ , i/?) = «a<i>, , H2 = {H+ , Hf^) = ej„vt^a ■ (1.3) 

This basis choice is known as the Higgs basis. In this basis, the 2HDM scalar potential becomes 
V = YiHIHi + Y2HIH2 + [Y3HIH2 + h.c] 

+ \Z^{H\Hif + \Z2{HlH2f + Z^{h\h^)(hIH2) + Zi{H\H2){HlH^) 

+ [\Z5{H\H2f + [Zq{h\Hi) + ZriHlH2)]HlH2 + h.c.} . (1.4) 



From the scalar potential in eq. (1.4) one can extract the squared-mass matrix for the three neutral states, 

/ Zi Re{Ze) -lm{Ze) \ 

M=v^\Re{Ze) i [Z3 + + Re(Z5)] + ^2/^^' -^MZg) . (1.5) 

yiMZe) -5lm(Z5) i [Z3 + Zi- Re(Z5)] + Y2/v^ J 

The field H2 and the coefficients Y3 and ^5,6,7 depend on the choice of basis, and thus are not physically 
meaningful by the mselves unless a symmetry (such as SUSY) is selecting a preferred basis. Thus, the mass 
matrix in eq. (1.51 is not invariant. This matrix is diagonalized with the transformation 

RMR^ = Md^ diag {mj , , ml) , (1.6) 
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with R defined as follows: "'^ 



ci2 -si2 0\ /ci3 -siaX no 

R = R12R13R23 = I Sl2 C12 1 C23 -S23 

1/ \s13O C13/ \0S23 C23, 

C13C12 -C23S12 - C12S13S23 -C12C23S13 + S12S23 ' 

C13S12 C12C23 - S12S13S23 -C23S12S13 - C12S23 
Sl3 C13S23 C13C23 



(1.7) 



where Cij = cos 0ij and Sij = aindij. The angles 9i2 and 6*13 are invariant; the third angle 6*23 depends on 
the basis choice. However, the product e^^^^H2 is invariant. Thus, one can express the squared-mass matrix 
entirely in terms of invariants: 

_ ( Zi Rc(Z6e-'«23) -Im(Z6e-*«^^) \ 

M = R23MR!^^=v^\ Rc(Z6e-'''23) YiciZ^e-'^'"^^) + /v"^ -\luv{Z^e-'^''^^^) , (1.8) 
V -Im(Z6e-'^^3) -ilni(Z5e-2'^^3) A^/v'^ j 

where is defined by 

A^ = Y2 + \\Z^ + Z4 - Re(Z5e-2»''^-')]w2 . (1.9) 

In this paper we use a mass ordering of the neutral fields such that m\ < 7712 < m^. The remaining Higgs 
particle has the mass 



mH± - ^JY2 + ^Z.v^. (1.10) 

In calculating the oblique parameters, certain invariant combinations of the mixing angles appear. These 
combinations are presented in Table |l] 



TABLE I: The U(2)-invariant quantities qkt are functions of the the neutral Higgs mixing angles 612 and 613, where 
Cii = cos Si,- and Si, = sin6'i, . 
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2. THEORY: THE OBLIQUE PARAMETERS S, T, AND U AND THEIR HIGHER-ORDER 

EXTENSIONS V, W AND X 

The oblique parameters S, T, and U provide an indirect probe of physics beyond the SM for theories 
with SU{2) X U{\) gauge content. They quantify deviations from the Standard Model in terms of radiative 
corrections to the gauge-boson two point functions, making use of the precise measurements available for 



^ In the CP-conserving 2HDM, only a single mixing angle is required. In this paper, the formalism of the CP-violating 
theory will be used, but as we will show, the numerical results for the oblique parameters do not change significantly in the 
CP-conserving limit. 
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parameters associated with W and Z boson resonances. New physics contributions is encoded in SHabiq"^), 
where a,b = j, , Z and 



(2.1) 



The 2HDM is a good candidate for analysis in the oblique correction formalism because the theory has no 
new electroweak gauge boson content and its couplings to light fermions enter only at two loops and are thus 
suppressed relative to the gauge boson couplings. 

If the new physics enters at the TeV scale, the effect of the theory will be well-described by expansion 
to linear order in q^, requiring only the three parameters (S, T, and U) originally defined by Peskin and 
Takeuchi g]: 



4s2 ' 



SUzzimD-SUzziO) 



<5n^^(o)-m;^(o) 



U 



4s2 
a 



snzz{ml)-snzz{o) 



(2.2) 
(2.3) 

(2.4) 



where Il'{q^) 



i-n((7^). These expressions, reprinted here from Burgess et al. [5], do not assume that 



the corrections vanish beyond linear order. If the vacuum polarization functions Ilaf, are expanded beyond 
linear order, six independent parameters are required to encode all the radiative corrections. The three new 
parameters may be defined as follows [5]: 



V = 
w = 

X = 



SU'zzir 



6nzz{ml)-6Ilzz{0) 



Sllwwim'^) - SllwwiO) 



( ^10 ) 



SUz^iml) ~ Snz.,{0) 



SU'z.iO) 



(2.5) 
(2.6) 
(2.7) 



These functions will be negligible if the new physics enters at a scale much larger than mz- 



2.1. Calculating the Oblique Parameters in the 2HDM 



Working on the assumption that the mass scale of the 2HDM is larger than niz but not so high that V, 
W, and X are irrelevant, we present here the contributions to the extended oblique parameters in the basis- 
independent formalism. These results are based on the one- loop calculations in rcfs. i6J and [2], with the 
modification that contributions from the A^^ A^H'^ and Af^Z^H^H" vertices have been added. The one- 
loop diagrams constructed from these two vertices cancel in the oblique parameters, but they are included 
here to show explicitly that (511^^(0) and (5n^-y(0) vanish as required [5]. The required vacuum polarization 
functions are as follows: 
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Slizziq^) = T 2 — 2 ( [^ii^o {q^,m%, ml) + q^^Ba {q^,m\,ml) + qI^Bq (g^, m|, ml)] 



, ml) 



+qliB22 (g^, m\,ml) + 921^22 {q^,m%,ml) + qliB22 {q^,m\,ml) + qliB22 {q^,ml, ; 
+qliB22 {q^,ml,ml) +931^22 {q^ ,m\,ml) + C2w^ B22 (g^ r7i^± , m^±) - ^[Aq {m\) + Ao {ml) 
+Ao (m^) +c2»2^o (m^±)] +to|So (g^^ml,™^) _ 5^2 (g^m|,m2)^^ (2.8) 



+911-^22 (g^, w|, TOi) + 921^22 (9^: "^L ™2) + 931^22 (9^) , mj) + qliB'22 {q^,'mj,ml) 
+qiiB22 {q'^,ml,ml) +qliB'22 {q^,m\,ml) +C2w^B'22 {q^ ,m]j±,m]j±) +mlB'g {q'^,ml,ml) 

-B'22{q\ml,ml)y (2.9) 



5^'wwi.q^) = [qiiB'o {q^,m\r-,m\) + qliB'„ {q^,m'^,ml) + qliB'^ (g^, m^, TO3)] 

+qiiB'22 {q'^,rn'^,ml) +(721^22 {q'^ , m'^ , mj) + qliB'22 {q^,m^,mf) +912^22 (q^m^±,ml) 
+922^22 {q^,fnjj±,ml) + ^32^22 {q^,m]j±,ml) + m^B'^ {q^,m'^,ml) ~ B^^ (q'^ , m'^ , mf)^ , (2.10) 



Snwwiq'^) = (^-"^^ [911 -Bo {q^,m^,ml) + qj^BQ (g^, m^, mg) + ggiBo {q^^rri^-, ml)] 

+qliB22 {q^,m^,m\) + qj^B 22 {q^,m'^,ml) +931^22 {q^,m'^,ml) + ^12^22 {q^,m'jj±,ml) 
+922-822 {q^,m]j±,ml) + ^32^22 ((7^,m^±,m3) - ^ [Ao(mi) +^0(^2) +-4o(to3) + Ao(m^±)] 



AvBq {q^,m^,m\) - B22 [q^ ,m^,m\)y (2-11) 

5n^^(g2) = ^ [B22 {q\mli,ml±) - \Ao (m^^i)] , (2.12) 

m;^(g2) ^ ^ [5^2 {q\m]j±,m]ji)] , (2.13) 

5nz,(g^) = [B22 {q\m\,,m\,) ~ \A, {m\^)] . (2.14) 



The following functions vanish, as required by the Ward identities, and thus may be removed from eq. (2.7) 
without loss of generality: 

STlz^fiO) = {B22 [0,TO^±,m^±]) - lAo (m^±) = 0, 

AUSyj Cyj 



m^^(O) = - [B22 (0,m^±,m^±)] - \Ao (m^±) = 0. (2.15) 
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3. WHAT THE 2HDM PREDICTS FOR THE OBLIQUE PARAMETERS 

3.1. Comparing with Experimental Limits 

In this section we compare the theoretical results for the oblique parameters with those of experiment. 
The experimental limits for S, T, and U are given in ref. [7]: 

5 = .014 ±.10, (3.1) 
r=.03±.ll, (3.2) 
;7 = .06±.10. (3.3) 

The limits for V, W, and X are known for a reference mass of 100 GeV ref. |H]: 

W^ = .ll±4.7, (3.4) 
y = .30±.38, (3.5) 
X = .38±.59. (3.6) 

To compare the experimental limits with the values predicted by the 2HDM, we randomly sampled the pa- 
rameter space of the scalar sector, choosing arbitrary values of^ Zi, Z3, Z3+Z4, Re{Zr, e~^'^^^), Re(Z6e~'^^^), 
lm.{Z^ g-2ie23^ g^]2(j l-[n{Zee~^^^^) , within the unitarity limits calculated in ref. [2]. Each set of parameters 
[Zi, etc.] determines the scalar masses and the values of the qke functions, producing unique values of S, 



T, U, V, W, and X via the expressions in Section 



2.1 



Since I2 is not constrained by unitarity, I2 



was used for all points. Although there is not a definite lower bound from experiment for mi as there is 
for the Standard Model Higgs mass, our code discards points that generate scalar masses less than m^.^ 
This cut-off was chosen so that our code would produce sensible results regardless of whether the oblique 
parameters are defined in the linear/quadratic approximation or defined to all orders in the expansion. 

The values of the oblique parameters generated in this random sample of 2HDM parameter space are 
shown in histogram form in Fig. |3.1[ With the exception of the T parameter, they fall well within the 
experimental bounds. Because the 2HDM generically produces values of T that are unrealistically large in 
magnitude, much of the parameter space will be eliminated on experimental grounds. This result will be 
explored further in subsequent sections. On the other hand, the distributions for S, U, V, W, and X fall well 
within the experimental limits, and thus are not useful in either constraining the parameter space or ruling 
out the 2HDM. Furthermore, these values all fall within roughly .1 of zero, implying that the 2HDM may 
be indistinguishable from the SM on the basis of these parameters. 

From the histograms, it is evident that the 2HDM predicts positive values for S and X. Since a random 
scan of parameter space tends to produce masses well above mz, the values for V, W, and X tend to be at 
least an order of magnitude smaller than for S, T, and U (more than one order of magnitude in the cases of 
W and X). This suggests that the 2HDM corrections are well described by S, T, and U; one does not gain 
much by going beyond the linear expansion approximation. These results are not specific to the particular 
value of I2 chosen; the allowed ranges of the Zi parameters (determined by unitarity) have more effect than 
the specific choice of I2 on the oblique parameters as long as I2 is less than or on the order of the square of 
the electroweak scale. Thus, choosing a smaller value for I2 does not produce larger ranges of the oblique 
parameters. 

The histogram for the X parameter has a different form from the others and merits some comment. In 
the 2HDM, X is a function of only, as illustrated in Fig. [2J To linear order in m|/m^± it is equal to 



^ We use Z-j + Z4 a,s a parameter rather than Z4 since Z4 is not by itself constrained by unitarity. 

^ Because the 2HDM produces very few points at low mass scales, discarding these points has an insignificant effect on the 
results. This subject will be discussed in more detail in SectionlSl 




^liiflf^iTTfh^WT^^ 



llllll^^llilllllliEi^^ 



200 
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50 



lllliinm. 



,v 



0.00002 0.00004 0.00006 0.00008 0.00010 O.OOOl. 



FIG. 1: Histograms showing the distribution of S, T, U, V, W and X values generated by a random sampling of 
2HDM parameter space. The experimental bounds are shown as hght gray shading {2a limits) and dark gray shading 
{la limits). 



^-1^ 60 (^-^^ 

Since X is proportional to l/m^j-, its lower bound depends on the maximum value of m/f±, which is 
controlled by and Y2. This is reflected in the histogram; the location of the sharp cut-off derives from 
the constraints from unitarity on the parameter Z^; if the unitarity bounds are relaxed, or if a higher value 
of I2 is chosen, the upper limit on m^fi increases, shifting the histogram data to the left. 

Other than X, the parameters fall into roughly bell-shaped distributions around their most probable values, 
exhibiting longer tails on the positive side than on the negative side. The statistics for these parameters are 
summarized in Table [TTl 

Because the experimental limits for each oblique parameter are not independent from one other, it is 
conventional to display results using plots of allowed ellipses in S-T space (or other 2D parameter spaces). 
The OPUCEM library [9^ can be used for this task. We adapted the OPUCEM algorithms for use in 
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FIG. 2: To first order, the parameter X is proportional to m^/m^±. 





Central Value 


Range 


Standard Deviation 


s 


0.030 


(-0.016, 0.19) 


0.029 


T 


0.19 


(-1.2, 8.5) 


1.6 


U 


0.00015 


(-0.010, 0.060) 


0.0088 


V 


0.00011 


(-0.0010, 0.014) 


0.0013 


W 


0.000033 


(-0.00033, 0.0086) 


0.00080 


X 


0.000016 


(8.3- 10-^ 0.00073) 


0.000066 



TABLE II: Predictions of the CP-violating 2HDM: Statistics for 1000 randomly produced points, with Y2 — m^. 



Mathematica to produce ellipses superposed on our scatterplots of data.''' The result for the S-T ellipse is 
shown in Fig.[3) ellipses in S-U and T-U space are shown in Fig.|4] These plots show that the points generated 
by the 2HDM are consistent with the experimental limits (with the exception of T, as noted previously). 




FIG. 3: Scatterplots of randomly generated points superposed on experimental limits in S-T space. The light gray 
shading shows the 1 a bounds and the dark gray shading shows the 2 a bounds. The second figure zooms in on the 
allowed region. 

The points with larger values of U that are evident in the S-U ellipse are not evident in the T-U ellipse 



* Our Mathematica program for generating the elliptical contours can be downloaded from the following site: 
|http://people. bridgewat er . edu/ ~ doneil/ ^ 
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FIG. 4: Scatterplots of randomly generated points superposed on experimental limits in S-U space [left] and T-U 
space [right]. The light gray shading shows the 1 a bounds and the dark gray shading shows the 2 a bounds. 



because there is a strong correlation in the 2HDM between U and T. Values of U that are larger than .01 
correspond to T > 1, and thus do not appear in Fig. |4] Incidentally, there is another strong correlation that 
crops up in the 2HDM: a near-linear relationship between V and W in the first quadrant. These correlations 
are illustrated in Fig. [Sj 




FIG. 5; Correlations appear in two sets of the oblique parameter values, T-U and V-W. The rest of the oblique 
parameters (not shown) do not appear to exhibit strong correlations with each other. 



These two relationships probably have a similar origin. For example, examining the definitions in eqs. (2.5) 
and (2.6), one notes that V (which is based on ZZ diagrams) is defined analogously to W (based on WW 
diagrams), so the fact that the vacuum polarization functions li^wio^) ^^nd TlzziQ^) have many similar 
terms leads V and W to have similar forms. To a lesser extent, the same phenomenon occurs with T and U. 
As a result of these correlations, the actual constraints on the 2HDM are more strict than those implied by 
the individual distributions illustrated in Figs. |3.1| 
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4. OTHER DEFINITIONS OF THE OBLIQUE PARAMETERS AND THE EFFECT OF THE 

g^-EXPANSION APPROXIMATION 



Two distinct definitions of the oblique parameters appear in the hterature; those of Peskin and Takeuchi [4), 
and the similar but physicially inequivalent expressions of Marciano and Rosner [10] . For the purpose of 
numerically analyzing the oblique parameters in the 2HDM, we have used in this paper the definitions of S, 
T, and U of Burgess et al. [5], which are generalizations of the MR definitions. They do not assume that 
the vacuum polarization functions Il{q^) are well-described by an expansion linear in . Kundu and Roy [8] 
advocate alternative expressions for S and U that also do not rely on the linear expansion approximation: 



S = 



U 



a 



snzzimD-snzzio) 



SUz^jml) 
'■z 



rar 



2 

■z 



5n 



WW 



(0) 



•'W 



6Uzziml)-SIlzz{0) 



SILzjiml) 



- s2 










) 









(4.1) 



(4.2) 



Since these expre ssio ns are based on the PT definitions, they have a slightly different physical interpretation 
than those of eqs. (2.2) and (2.4 1. To linear order in q^, there is no difference between the two, but beyond 



linear order they differ by a factor on the order of the X parameter [8] . To assure the reader that no significant 
errors are introduced by applying the PDG experimental limits (based on the PT definitions) to theoretical 
computations in the Burgess scheme, comparisons of the two sets of oblique parameters are shown in Fig. |6] 
To examine the behavior of the functions S and U as the mass scale of the new physics (NP) increases, we 
define 



M 



NP 



\ (toi +7712+7713+ mii± ) , 



and choose 



Re(Z5e-2«^23) 



Zi = .5 (Z3) + .5, 
: Re(Z6e-'''=^3) = -.1, 



5^2 = mij^, Z4 = 0, 

Im(Z5e-2«^23) ^ lm{Zee- 



= .1. 



(4.3) 



(4.4) 



With these values for the scalar couplings, dialing up Z3 simultaneously increases the neutral masses 777.^ 
and the charged Higgs mass 777 . Thus, one can analyze the dependence of the oblique parameters on the 
mass scale M^rp by varying Z3. The results in Fig. [6] show that the two definitions converge in the limit 
q^ <tC M^p, as they should. Even at lower mass scales the difference is insignificant for S. There is noticeable 
difference at low Mj^p for U, but the difference is dwarfed by the experimental error in U [see eq. (3.3)], so 
this should not undermine confidence in our numerical results. 

Meanwhile, there are also different ways of calculating W and V. The definitions advocated by Kundu and 
Roy in ref. [8] are mathematicall y eq uivalent to those introduced by Burgess et al. [see eqs. (2.5) and (2.6)], 
so the experimental limits in eq. ( |3.6| ) apply equally well to either. For completeness, one can also check that 
these expressions for W and V converge at high mass scales to those defined in the quadratic approximation: 



n(g2) = n(o) + q'^u'io) + h^u"{o). 



In this approximation, V and W become: 



w 



6nzz{ml)~SIlzziO) 



^2 



SWww{m\f) - OTivM/(0) 



777 



^^wwi^) 



w 



(4.5) 



(4.6) 



(4.7) 



In reproducing these equations from ref. [5j we have converted them to our notation, in which n'(q^) = 
-£pH{q^). Figure jrj show that V and W in this quadratic approximation converge to the non-approximate 
forms as the mass scale of the new physics increases, as expected. 
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FIG. 7: Comparison of V and W in the formulation by Burgess et al. and in the quadratic approximation. 



5. CONSTRAINING THE 2HDM WITH THE OBLIQUE PARAMETERS AND UNITARITY 

BOUNDS 



5.1. Constraints on the Mass of the Charged Higgs Boson 



In this section we discuss what the obhque parameters indicate about the charged Higgs mass mf^± . It 
fohows from eq. (1.10 1 that the unitarity bound \Z^\ < 8n affects mu±. In Fig. [Sj the allowed region for 
rriHi is shown. Without knowing Y2, one cannot put a definite upper bound on m^-i. However, by using 
the experimental limits on the oblique parameters, a more precise bound may be obtained. 




200 400 



500 SOO 1000 



4y; 



FIG. 8: Range of rnu± values allowed by the unitarity limits on Z3, as a function of 
^/Y2 around the electroweak scale, the upper limit on m^± is about 900 GeV. 



/Y2 (in GeV). For values of 
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The most promising of the obhque parameters in this regard are S and T. For the purposes of this analysis, 
we will use the PDG values of S and T with U fixed to be zero, since most values of U were found to be 
within .01 of zero in section ISTTl 

S = 0.03 ±0.09, 

T = 0.07 ±0.08. (5.1) 
The behavior of S and T as a function of m^± is shown in Fig. |9] The points in the scatterplots were 



generated as in section 3.1 with the added inclusion of extra points (shown in gray) that appear when the 
unitarity constraint on is relaxed by a factor 4 (ie, jZgl < 327r instead of jZaj < Sir). The plots show 
that low values of S are associated with higher values of ■ The distribution of S values predicted by the 
2HDM (according to randomly selecting the scalar parameters within the unitarity limits) has a peak around 
S = .02 (a value that agrees well with experiment), which correlates to a relatively high value of charged 
Higgs mass {mH± > 600 GeV) . The results for T also appear to disfavor low values of mH± ; although points 
producing m//± < 250 GeV are easily generated, they tend to produce unrealistically high values of T. It 
appears, then, that the 2HDM predicts 

mH± Z 250 GeV, (5.2) 

with higher values favored. 




FIG. 9: Scatterplots of randomly generated S and T values, with Y2 = m^. The black points are allowed; the gray 
points violate the unitarity limits on Z3. The experimental values are shown as a dashed line; the gray shading shows 
the 1(7 and 2a error bars. 



One should be careful in interpreting the results in Fig. |9] For one thing, the experimental limits in 



eq. (5.11 assume a Standard Model Higgs mass of = 117 GeV. In our analysis, we have been fixing 
TO0 = mi. Since mi does not equal 117 GeV in general in this analysis, the experimental bounds will 
not be completely accurate for all points. To determine whether this is significant, let us examine the S 
and T scatterplots in more detail. Essentially what appear in the scatterplots are sets of possible curves, 
which exceed the experimental bounds at different points. One such curve is shown in Fig. |5.H generated 
by choosing 

Zi=3.5, ¥2 = ml., ^3 + 2-4 = 4, 
Im(Z5 e-2*^^-^) = Re(Z5 £-2*^23) ^ ReiZec-'"^-') = lm{Zee-''^^^) = .5, (5.3) 

with varying values of Z3 to produce a function of mjj± . 

The neutral Higgs boson masses for this choice of parameters are 

mi = 321 GeV, ma = 363 GeV, mg = 483 GeV. (5.4) 

For rricf, = 300 GeV, the PDG gives different experimental limits: 

S = -0.07 ±0.09, 

T = 0.09 ±0.08, (5.5) 
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FIG. 10: The behavior of S and T at the point given by eq. ( |5.3[ ), where mi — 321 GeV. With this value as a reference 
mass, the experimental limits shift slightly [see eq. (5.5|]. 



SO the limits on the plots have been adjusted accordingly [see Fig. 5.1 . Thus, it does not appear that shifting 
the experimental limits to correspond to a higher SM reference mass appreciably affects the conclusion of 
eq. (IK2I). 



However, one more caveat is required, which we will illustrate as follows: Consider the following parameter 
values. 



Zi = 3.5, Y2 = m'^, ^3 + ^4 = 1, 



(5.6) 



Im(Z5 e-2*«=3) = Re(Z5 e~^''^'^) = Re{Zee~''^^^) = lin{Zee-''^^-') = .5, , 
which generate the following mass values: 

TOi = 117 GeV, ma = 221 GeV, mg = 472 GeV. (5.7) 
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FIG. 11: The behavior of S and T at the point given by eq. (5.6 1, where mi — 117 GeV, 



The behavior of S and T at this point in parameter space is shown in Fig. 11 Note that in this figure the 
T curve has shifted to the left, into a region not represented on the scatterplots of Fig.[9j In other words, the 
lower bound in eq. (5.2) may be too restrictive, since it represents mostly higher sets of neutral Higgs masses. 



Thus, while one can conclude that the 2HDM generically predicts values of mff± that are above 250 GeV, 
additional information about the theory may render the constraints from Fig. |9] inapplicable. (An example 
would be if the lightest Higgs mass was detected below 150 GeV.) It should be noted, however, that the point 
defined in eq. (5.6) is particularly nongeneric; it was tuned to give an especially low mi value. In general, the 
2HDM predicts mass scales that are fairly high (400 - 800 GeV), as is evident from the distributions shown 
in Fig 12 To demonstrate that this favoring of high scales is not an artifact of our choice to discard points 
with mi < mz, the actual distribution produced by our program is compared in the figure to one with all 
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real- valued masses allowed. The figure shows that imposing this minimum mass condition introduces only a 
very slight biasing in the high-mass direction. 
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FIG. 12: Mass scales that result from a random sampling of 2HDM parameter space at Y2 — m^. The mass scale 
Mjvp is the average of the three neutral masses and m^^i [see eq. (4.3 1]. The cut-off at the high end is due to 
imposing unitarity limits. The strong overlap between the distributions shows that the lower limit is not due to 
imposing rrii > mz (as we have done throughout this analysis). 



A final comment on the charged Higgs mass: If the parameter X was known more precisely from experiment, 
it would be possible to determine from X alone. X, as can be seen from eq. (3.7), is an analytic function 
of mj^± and does not depend on any other scalar masses or couplings. Unfortunately, the experimental 
bounds are 2 orders of magnitude too large to be useful in determining m/^± (see Fig. ^ . 



5.2. Constraints on the Splitting between the Charged and Neutral Scalar Masses 



The difference between the charged Higgs mass {rau±) and the heaviest neutral mass (rris) in the 2HDM 
is strongly constrained by the experimental bounds on the T parameter. In Fig. |5.3[ the behavior of T at 
the point defined in eq. (5.3) is shown. The experimental bounds on T constrain to be higher than mfj±. 
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FIG. 13: The T parameter at the point given by at the point given by eq. ( |5.3[ ), with mi = 321 GeV. The gray 
shading shows the 1 a and 2 a error bars. The experimental limits are shifted to match a 300 GeV reference mass, 
as in eq. (5.5). 



that for Yy 



the 



2HDM seems to strongly favor mu± < m^. In particular, the only points within the 2cr bounds on T are 



In fact, this behavior is fairly generic-one can see from the scatterplots in Fig. [T4| t 

n tb 

those with 77i^± — 7713 between —600 and +50 GeV. 

The role of T in constraining the splitting becomes especially important for higher values of ¥2- For 
Y2 — (ITeV)^, there are plenty of points for which TO//± > ma, but only those for which m//± is less than 
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FIG. 14; Scatterplots of randomly generated T values, for Y2 = (mwY P^ft] and Y2 — (ITeV)^ [right]. The gray 
shading shows the 1 a and 2 a error bars. 



100 GeV greater than 7713 survive the experimental constraints on T (see righthand image in Fig 14). Thus, 
allowing for a wide variation in the possible values of Y2 , one can conclude from this analysis that 



600 GeV < mH± - mg < 100 GcV. 



(5.8) 



Again, one should check that this conclusion remains valid for a light value of mi. Returning to the special 
point defined in eq. (5.6), with mi = 117 GeV, one sees from Fig. 15 that the curve has shifted entirely into 



the m3 > mH± region. This is consistent with the result in eq. (5.8) 




FIG. 15: The behavior of T at the point given by eq. (5.61, with mi = 117 GeV, 



6. THE CPC MODEL VS. THE CPV MODEL: DO THEY MAKE DIFFERENT PREDICTIONS 

FOR THE OBLIQUE PARAMETERS? 



The preceding analysis was done in the CP- Violating 2HDM. To see if the results would change if the scalar 
sector of the model conserved CP, the oblique parameters were calculated in the CP-conserving limit, taking 
Im(Z5 = Im(Z6e~*^^^) — 0. This produced small shifts in the distribution of the oblique parameters; 

two examples of which are shown in Fig. 16 The differences between the CP-violating and CP-conserving 



results were in all cases much smaller than the width of the distributions, confirming that the numerical 
findings in this paper apply to the CP-conserving 2HDM as well. Unfortunately, this result precludes the 
possibility of predicting from experimental values of the oblique parameters whether the 2HDM would have 
CP-violation in the scalar sector. 



CONCLUSION 



This numerical analysis of the CP-violating 2HDM shows that the model generically predicts large scalar 
masses, with average masses in the 400-600 GeV range. It is unsurprising, then, that a random sampling of 
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FIG. 16: Taking the CP-conserving limit shifts the distribution of S values slightly in the negative direction [left] and 
shifts the distribution of T values slightly in the positive direction [right]. The experimental value of S and T are 
shown as dashed lines; the light gray shading shows the 1 a error bars. 

the parameter space produces values of the extended oblique parameters W, V, and X that are at least an 
order of magnitude smaller than those of the traditional oblique parameters S, T, and U, since the extended 
oblique parameters are expected to be small except where the mass scale of the new physics is very close to 
the Z boson mass. Larger values of W, V, and X are consistent with the 2HDM, but are not characteristic 
of the model. We also find that the 2HDM produces values of U within .01 of zero, suggesting that the 
behavior of S and T is best compared to the experimental values where U has been fixed to be zero. The 
T parameter, for which the 2HDM produces a wide range of values, puts constraints on the possible values 
of the charged Higgs mass and even more so on the splitting between the heaviest neutral scalar mass and 
the charged Higgs mass. The 2HDM favors large mu± values, of at least 250 GeV, and even larger values 
of TO3, the heaviest neutral Higgs boson. In particular, we find —600 GeV < m^i — ma < 100 GeV. This 
suggests that if a light scalar particle (say ^ 120 GeV) is discovered at the LHC, the possibihty remains of 
heavier scalars being found closer to the TeV scale, with the charged Higgs particle expected to be lower in 
mass than the heaviest neutral state. 

These numerical results are not significantly affected by altering the definitions used for the oblique 
parameters, nor by going to the CP-conserving limit. 
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